Impossibility of large phase shifts via the "giant Kerr effect" with single-photon 

wavepackets 
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An approximate analytical solution is presented, along with numerical calculations, for a system 
of two single-photon wavepackets interacting via an ideal, localized Kerr medium. It is shown that, 
because of spontaneous emission into the initially unoccupied temporal modes, the cross-phase 
modulation in the Schrodinger picture is very small as long as the spectral width of the single- 
photon pulses is well within the medium's bandwidth. In this limit, the Hamiltonian used can be 
derived from the "giant Kerr effect" for a four-level atom, under conditions of electromagnetically- 
induced transparency; it is shown explicitly that the linear absorption in this system increases as the 
pulse's spectral width approaches the medium's transparency bandwidth, and hence, as long as the 
absorption probability remains small, the maximum cross-phase modulation is limited to essentially 
useless values. These results are in agreement with the general, causality- and unitarity-based 
arguments of Shapiro and co-workers. 



I. INTRODUCTION 

It was suggested by Chuang and Yamamoto [l[ that the 
cross-Kerr effect (or cross-phase modulatation) between 
two optical fields could be used for so-called "dual-rail" 
quantum logic with photonic qubits, provided sufficiently 
large nonlinearities could be generated. The goal is to 
be able to perform a transformation like the following 
(where and 1 refer to the number of photons in the two 
interacting modes): 

-L(|0) + |l))|l)^i=(|0) + e -^l))|l) (1) 

For cross-phase modulation by a x^ medium, the phase 
shift <f> = nn\ri2, where n\ and ni are the number of 
photons in the respective modes, and k is a constant. If 
k can be made as large as tt, then Eq. JT]) shows that we 
have a conditional logical operation which, in the basis 
of states (|0) ± |1))/V2, is equivalent to the CNOT gate. 

As already pointed out in [2], however, the very large 
Kerr nonlinearites needed to achieve this sort of phase 
modulation are, in ordinary materials, always associated 
with large absorption losses. Over the years, a number of 
methods have been suggested to overcome this difficulty, 
mostly revolving around the use of electromagnetically- 
induced transparency, or EIT, to eliminate the linear 
absorption and, at the same time, increase the nonlin- 
ear dispersion of the medium 0, 0, @|. Nevertheless, in 
2006, J. H. Shapiro published a study [1] showing that 
the causal, non-instantaneous behavior of any x^ non- 
linearity would always prevent the transformation ([I} 
from happening with high fidelity for large <fr. Cen- 
tral to Shapiro's study was a multimode treatment of 
the two propagating one-photon wavepackets, something 
that most previous studies had not considered. 

Shapiro's original argument was very general, and in a 
later publication Q it was stated that it was not imme- 
diately apparent whether it applied to the EIT schemes. 
One such scheme that had, in fact, received a multi- 
mode quantum-field treatment was that of Lukin and 



Imamoglu @ (see also the discussion in Q), based on 
the so-called "giant Kerr effect," originally proposed by 
Schmidt and Imamoglu 03 (which is, in one form or 
another, at the heart of all the EIT proposals). The 
conclusions of 8, 9] are, however, somewhat ambiguous, 
because they suggest that very large phase shifts are pos- 
sible, but at the cost of large changes to the modal struc- 
ture of the pulses. 

The goal of the present paper is to show that with an 
appropriate, idealized, but local Hamiltonian, that repro- 
duces the Heisenberg-picture evolution equations of refs. 
d, Q , large phase shifts in the Schrodinger picture (as in 
Eq. {TJ)) are, in fact, impossible, if the traveling fields are 
described by localized, single-photon pulses. This is in 
agreement with Shapiro's prediction for the "fast nonlin- 
earity" case. In fact, the analysis presented here clearly 
shows that a large phase shift is only possible in the limit 
in which the pulse's spectral width approaches the band- 
width of the nonlinear medium. In that case, however, 
it is also shown here, for the specific atomic configura- 
tion considered [Io| in the derivation of the giant Kerr 
effect, that EIT becomes ineffective and the losses due to 
absorption in the medium approach unity. 

The paper is organized as follows. The instantaneous 
response limit, with an idealized Hamiltonian, is consid- 
ered in Section 2. The actual time response of an EIT 
medium is considered in Section 3, where a relationship 
between absorption losses and the pulse's spectral width 
is derived. Section 4 contains some further discussion 
and conclusions. 



II. INSTANTANEOUS RESPONSE LIMIT 



A. Hamiltonian, locality, and Heisenberg picture 
results 



Suppose that, somehow, one has managed to produce a 
medium that leads, to a sufficiently good approximation, 
to a field evolution described by the following Hamilto- 
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Here i? Q and E\, are Schrodinger-picture field operators 
given by 
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with i5a "* their Hermitian conjugates. The quan- 

tization volume has cross-sectional area A and length L. 
The total number of modes to be added, for each field, is 
M = 2n max + 1, and it is determined by L and the band- 
width of the optical medium, Aw me dium = 2ttcM / L. The 
number of initially-occupied field modes may be much 
less than M, as will be discussed below. The following 
Section discusses how the "giant Kerr effect may 
lead to the kind of interaction Hamiltonian that appears 
in ^ under appropriate conditions. 

An important property of the Hamiltonian ^ is that 
it is local. Suppose that the field state is described by a 
pair of wavepackets, which for simplicity we will take to 
be identical: 
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Here |l„) a is a state with one photon in the n-th "a" 
mode, and zero photons in all the other modes, and sim- 
ilarly \l n )b- The field intensity in such a state is given 
by 



I a {z) = (E^E^) 
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and the expectation value of the energy is 
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n\c n \ 2 + hel I a (z)I b (z)dz (6) 



For pulses with a symmetric spectrum, such as those 
which will be considered here, the first term in ([6]) van- 
ishes. The second term, on the other hand, which rep- 
resents the interaction energy between the field and the 
material medium, vanishes if (and only if) the pulses are 
not inside the medium. This is a physically reasonable 
requirement for any Hamiltonian that one might want to 
use to describe the interaction of a finite pulse with a 
localized medium. 



If the Heisenberg equation of motion is used with 
the Hamiltonian @, one gets the following propagation 

equations for the field operators E^ (t, z) and E^'(t, z) , 
in the Heisenberg picture: 



dt dz 



ot az 



(7) 
(8) 



where k = e(hu)Q / 'cqA) , and the right-hand side is as 
shown for zq < z < zq + 1, and elsewhere. This follows 
from the assumption that the bandwidth of the medium 
is large enough to justify the approximation 
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Strictly speaking, this requires also that the bandwidth 
of the pulses be smaller than the bandwidth of the inter- 
action, an assumption to which we will return shortly. 

Eqs. ([HI) are exactly the same as obtained by Lukin and 
Imamoglu in [8| , minus all the extra complications arising 
from the different group velocities in the medium and 
the wavepacket compression. These complications will 
simply be ignored here in order to concentrate on the 
basic difficulties caused by locality and the multimode 
nature of the field. As explained in the Eqs. ([5]) 
can be self-consistently solved by integrating along the 
characteristics to get, at any point z within the medium 



,(+) , 



x exp 



i — iz 

c 



z )Eti(t>,zo)Ett(t>,zo) 



(10) 



with the local time t' = t — (z — z )/c. Eq. (JTUJ) 
can be verified by direct substitution in ((SJ), not- 
ing that it implies the equality E^Ttf (t, z)E^ b ' (t, z) = 

El~b(t', z o)E^(t', zq). For z > z a + I, one can use 
free propagation backwards, and the fact that the field 
only undergoes multiplication by a unitary operator, to 
rewrite the result in terms of the t — operators: 
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At first sight, Eq. (fTTj) might seem to be exactly what 
we want, since it suggests that each of the two fields 
acquires a phase that is proportional to the intensity of 
the other one. The fact that the actual phase shift appar- 
ently depends on the local intensity, at different points in 
the wavepacket, may be slightly worrisome, but Eq. 
at least suggests that nothing should prevent one from 
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making the phase at, say, the center of the wavepacket, 
as large as one might want to. The situation looks quite 
different, however, in the Schrodinger picture, to which 
we turn next. 



where, for convenience, the parameter rj = e(hujQ / €qAL) 2 
has been defined. One can next introduce a new set of 
coefficients, U/x(i), defined by 



B. Time evolution in the Schrodinger picture 

In the Schrodinger picture we write the state of the 
system as the double sum 

i^w>-EE c »( i ) r2l,< " +m)ci/L i 1 »>«i 1 ™> i ( 12 ) 



where the coefficient c nm (0) (two indices) equals the 
product c n (0)c m (0) (single index) at t = 0. The equation 
of motion for c nrrl is 



2n m „ x — M 



E 



(15) 



v— Irirn 



Note that in (fTS} . as in all other sums over v for constant 
H, the index v increases in steps of 2, so there are 2 ! Yiwiax 

+ I terms in the sum (and \i ranges from — 2n max 
to 2n max ). The obey the equation of motion (easily 
derived from (jT4j) ) 
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This can be integrated analytically under some approxi- 
mations that are equivalent to the ones used in the previ- 
ous section. To begin with, introduce a new set of indices, 
fi and u, that stand for the sum and difference, respec- 
tively, of n and m. Then c„ M = c nm with n — (jj, + v)j1 
and to = [p, — v)/2, and we have 
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Equation (|16|) can be integrated by introducing an enve- 
lope function f[t,z), defined as 



uj(ct-z) I L 
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(14) which satisfies 
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It is easy to see that, for large enough n max , the expres- the previous section, one finds, for z inside the medium, 
sion in square brackets in (fTS)) converges to (2n max + 
l)L5(z — z'), which means that we have 



' d_ 
dt 



f = -ir)MLf(t, z), z < z < z + l 

(19) 

where, again, M — 2n max + 1 is the total number of 
modes, and the right-hand side of (|19[) vanishes outside 
the medium. Integrating along the characteristics, as in 



f(t, z) = e -*J^£(»-*o)/c/(o, z - ct) 

= e - tv ML(z-z )/c^2 v ^Q) e ~^ 



ifi(ct-z) j L 



(20) 
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which can now be substituted into (fl4|) (via the defini- 
tions (JTSJ) and ([17])). to yield 

<w = - «?e v (o) 

A*' 

/■ZO+i 

x / g-STri^ -fi)(ct-z)/L e -ir]ML(z-z )/c ^ 

(21) 

The expression on the right-hand side of (f2Tjl can now 
be directly integrated. For the purpose of comparing the 
state of the field after the interaction to the state be- 
fore the interaction, it is convenient to concentrate on 
the value of the coefficients c„ M at the time T = L/c (i.e., 
the quantization time), at which point, in the absence of 
interaction, the pulse should return to its initial state, 
since the traveling- wave formalism we are using is equiv- 
alent to periodic boundary conditions. In that case, the 
integration of (|2"Tj) over time from t = to t = T selects 
only the ji = fi term in the sum, and we have 

c M „(T) = c^(0) + (e-^ MLl ' c - l) MO) (22) 



C. Fidelities, and numerical results 



Equation (|22|) can be used to calculate the overlap be- 
tween the initial state and the state at the time T, from 
which, in turn, a number of other useful results can be 
derived. Defining, for simplicity, 
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(where the last expression uses k as defined in the previ- 
ous section, for comparison with the Heisenberg-picture 
result, Eq. (fTTj) ). we have 



(V(o)|VCO) = i- 
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Recall that the original goal (Eq. |T])) was to leave the 
original two-photon state invariant except for a phase 
shift. The fidelity T§ is a measure of the success of this 
operation. Note that if \v^(0)\ 2 /M = 1, we have 
To = 1 an d (f> — 3>. It is important, therefore, to calculate 
this quantity. Note that the expression ([5]) for the single- 
wavepacket intensity can be rewritten as 



I a (z) = (EMEW) 
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e AL 



5>;(0)c m (0> 



from which it follows that 
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and therefore 
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I a (z)dz = - 



Lj L P a (z)dz 



So ' I a {z)dz 
(27) 

This can be related to the pulse bandwidth as follows. 
First, note that the assumption of a localized pulse means 
that it is legitimate to extend all the integral arguments 
in (|27[) from minus infinity to infinity. Then, introduce 
the spatial Fourier transform P(k) of the function I a (z), 
so that 
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Then we have 
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(28) 



(29) 



Now suppose that the function P(k) is peaked at k = 
(as it should be if we have correctly separated the slowly- 
varying part of the pulse form its carrier frequency) , and 
that it is negligible outside of an interval of width Afc. 
Then |P(fc)| 2 /|P(0)| 2 < 1 for all k, and therefore 



Em°)i : 



< —Afc 
~ 2tt 



(30) 



where the equality holds only for a "square" P(k) (con- 
stant in the interval Afc, and zero outside of it); but this 
corresponds to a spatially non-localized pulse, whose in- 
tensity decays only as 1/z. We can then assume that 
(|30| is always a strict inequality. Returning to our origi- 
nal formulation in terms of M discrete modes spaced, in 
frequency, by 2nc/L, we conclude that 
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pulse 
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(31) 



which should always be less than 1 (note that Aw pu i se has 
been defined through the effective support of P(k), and 
will typically be larger than the conventional "standard 
deviation" of lj for the wavepackets considered) . In terms 
of the quantity r, we can express the fidelity and the 
phase <j) as 



and 



tan 



r sin $ 



r) 



1 



r cos $ 



(32) 



(33) 



Note that we can always make the distortion of the orig- 
inal wavepacket negligible by letting r — > 0, but this is at 
the expense of an extremely small phase shift. The op- 
posite case, r — > 1, also leads to a high fidelity, this time 
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taken to start at zq = L/2 and had extension I = L/2. 
The agreement with the analytical approximation is bet- 
ter for the broader pulse because the replacement of the 
term in square brackets in p8[) by a delta function is 
a better approximation in that case. The calculations 
show that the theory may underestimate somewhat the 
achievable phase (j>, but large values of <fi still correspond 
to very small fidelities. 



D. Discussion 



FIG. 1: Parametric plot of the fidelity To versus the phase 
shift <$> for r = 0.4, 0.3, 0.2, as $ is varied from to 2ty. Dots: 
result of numerical calculations with a Gaussian pulse, with 17 
modes and a/L — 0.059 (corresponding to r = 0.4, outermost 
dots), and a/L = 0.078 (corresponding to r = 0.3, innermost 
dots); the numerical calculations cover only the range < 

<1> < 7T. 



with potentially large phase shifts, but it is forbidden by 
locality. By this it must be understood that it is simply 
not possible to make r arbitrarily close to 1 with a local- 
ized pulse. For instance, suppose that the intensity I a is 
proportional to a Gaussian e _< - z_Zl - 1 I" . The right-hand 
side of P?]) then evaluates to 1/ (er V27r) , and we have 



L/a 1 
M J 7 ^ 



< 0.4 (Gaussian) 



(34) 



since, in order to describe the pulse adequately, the 
number of modes M must be at least of the order of 
L/a (in other words, using fewer modes results in a 
non-localized pulse). Similarly, for a hyperbolic secant 
Ia(z) ~ sech(z/er), we find 



L/a 2 



< 0.2 (hyperbolic secant) 



(35) 



Since order of magnitude arguments are always uncertain 
regarding such things as factors of two, however, and the 
temptation to look for an "optimal pulse shape" is strong, 
it is important to keep in mind the absolute bound (f3"Tj) . 
and also the considerations to follow in the next section, 
which will show, for the specific example of a "giant Kerr 
effect" medium, how things deteriorate when one tries 
to fit the pulse's spectrum too tightly in the medium's 
transparency window. 

Parametric plots of To versus <j> for various values of 
r < 0.5 are shown in Figure 1. Note that, when r = 1/2, 
<p = <&/2 and To — cos 2 (<f>/2). The figure also shows 
the results of numerical calculations (based on the di- 
rect integration of Eq. (fT3"|) ) for M = 17 modes and two 



Gaussian pulses with I a 



-(z-^f/a* 



L/A, and 



a = 0.059L (corresponding to r — 0.4) and a — 0.078L 
(corresponding to r — 0.3). The nonlinear medium was 



The results in the previous subsection indicate that the 
parameter r that determines the maximum achievable 
phase shift decreases, as the ratio of the pulse's spectral 
width to the medium's bandwidth. 

One may wonder why, once enough modes have been 
included in the calculation to describe the wavepacket 
properly, the addition of more, empty, modes should 
have a degrading effect on the performance of the sys- 
tem. The answer lies in spontaneous emission. Consider- 
ing the action of the Hamiltonian Q on an initial state 
with only one photon in each pulse, the two sets of an- 
nihilation operators, and will first produce a 
uniform vacuum state, and then the creation operators 



E { - ] 



and Ei ' may replace the a and b photons into 
any of the available temporal modes, regardless of where 
they may have been originally. There are some momen- 
tum and energy conservation constraints, enforced by in- 
tegrals over position and time respectively, but as (for 
instance) Eq. (|14[) shows, an initially unoccupied pair of 
modes with an arbitrary v (difference between "a" and 
"b" temporal frequency) can be created out of any of 
the preexisting pairs of modes with the same /i (sum of 
the "a" and "b" frequencies), without incurring any en- 
ergy or momentum penalty. This is also apparent from 
Eq. (J221). 

This problem is particularly acute for single-photon 
wavepackets. For pulses containing an appreciable num- 
ber of photons, n, the action of the annihilation opera- 
tors results in modes that are still highly populated, and 
so stimulated emission will take place preferentially in 
those modes. In other words, one expects this problem 
to decrease as 1/n (the ratio of spontaneous to stimu- 
lated emission) , as the number of photons in the pulses 
increases. 

To get the single-photon case to work, one might con- 
template modifying the Hamiltonian, so that, for in- 
stance, instead of the negative-frequency field operators 



Ea and E^ ' one would have a weighted sum of creation 
operators, more closely matching the spectrum of the 
incoming pulse (of course, by Hermiticity, the positive- 
frequency field operators E^ and E^ would also have 
to be modified). This amounts to introducing some of 
the effects of dispersion in the medium, but it cannot be 
done arbitrarily, since there are physical rules (such as the 
Kramers-Kronig relations) that govern these things. In 
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FIG. 2: Level scheme for the giant Kerr effect. E a and E), 
are weak (in this paper, single-photon) fields; Q c is the EIT 
"coupling field." 



particular, strong dispersion is typically associated with 
absorption. As will be shown in the next section, even 
the extremely weak residual absorption present in the 
giant (EIT-enhanced) Kerr effect is enough to prevent 
one from taking the limit r — > 1 in the results presented 
above. 

It may be worth considering for a moment the extreme 
case of a "toy" Hamiltonian that would work with any 
pulse shape. This could be achieved by replacing the 
interaction part of (J2J> by 



Hi = file I aia n 



(36) 



Unlike |(3J|, this Hamiltonian does not create any pho- 
tons in initially unoccupied modes; yet, it is unphysi- 
cal, because it is completely nonlocal: the pulse will be 
interacting with the medium wherever it might happen 
to be. This highlights the importance of doing multi- 
mode quantized-field calculations properly. Formally ei- 
ther one of the Hamiltonians @ or (|3"6")l could be con- 
sidered as a possible generalization of the single-mode 
Kerr Hamiltonian a^ab^b, but they yield very different 
predictions in the multimode case, and only one of them 
is (approximately) physical. It is to the terms neglected 
in this approximation that we turn in the next Section. 



III. GIANT KERR EFFECT AND MEDIUM 
BANDWIDTH 

In the previous section it was shown that in order 
to achieve a relatively large phase shift one should try 
to make the pulse's bandwidth as close to that of the 
medium as possible. However, when one does that, the 
medium's absorption is not negligible anymore. 

A way to approximately realize the Hamiltonian @ 
is by making use of the "giant Kerr effect" introduced 



in [HI. The scheme (illustrated in Figure 2) makes use 
of electromagnetically-induced transparency, or EIT, to 
open a "transparency window" in the medium for the 
field E a (the field E does not need it, since the level 2 
will typically be unoccupied, and the detuning Af, will be 
taken to be large), as well as to enhance the dispersion 
of the medium and with it the Kerr nonlinearity. As 
shown in [ffj , if the optical Bloch equations for the atom 
are solved in steady-state, under the assumption that 
absorption is small, one obtains for the atomic dipolc 
amplitudes in the two transitions the result 



Pa = I 



Pb = i 



4df 3 d| 4 



^2( 74 / 2 + i A h ) 



\E h \ z E: 
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E„\ 2 El 



(37) 
(38) 



where G?i3 and <i 2 4 are the dipole matrix elements for 
the two transitions. Multiplying each dipole amplitude 
by the corresponding field, and adding the contributions 
of all the atoms by integrating over the spatial extent 
of the medium, then yields an interaction energy of the 
form ©, with 



_ A 



(39) 



where p is the volume density of atoms in the medium 
and A the cross-sectional area of the beam, and the as- 
sumption Ah 3> 74 has been made. As discussed in, 
e.g., Q, this steady-state approximation, appropriate 
for a continuous-wave field, neglects a number of impor- 
tant dispersive effects that result in a slowing down and 
broadening of the E a pulse. These complications were 
discussed by Harris and Hau in and possible ways 
around them were suggested by Lukin and Imamoglu in 
[1]. Here the propagation effects will be ignored, in or- 
der to concentrate, instead, on the consequences of the 
temporal variation of the pulse at the location of each 
atom. 

For an optically dense medium, the transparency win- 
dow is a Gaussian of width 



AtJtrans — 



o 2 



1 



V r 3i73i y/po-J 



(40) 



(see, e.g., [TT() where a a — 3A 2 /27r is the on-resonance 
absorption cross-section of the atom in the 1 — > 3 tran- 
sition; is the spontaneous emission decay rate from 
level 3 to level 1, and 731 > T3i is the total decay rate of 
the (1,3) coherences, including dephasing and decay to 
other levels (such as 2). The residual absorption inside 
this window has been discussed using a purely semiclas- 
sical treatment in [l3j]; here, for single-photon pulses, it 
will be estimated (and with it, implicitly, the width of the 
window itself) as follows. Considering only a one-photon 
wavepacket in field a, and ignoring field b altogether for 
simplicity, an initial state J2 n c n(0)\l n ) a\^-} (where the 
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second ket refers to the atomic state) can evolve into a 
superposition 



\m) 



c n (t)e- m " t \l n ) a \l)+C 2 (t)\0) a \2)+C 3 (t)\Q) a \3) 

(41) 



(here and in what follows, to = 2ttc/L), under the non- 
Hermitian Hamiltonian 



H appr = J2 hn ^ a - d i3 (E ( a +) 1 3) <1| + 4 _) 1 1) (3 1 



2 (|2)(3| + |3)(2|)-^|3)(3| 



(42) 



For simplicity (in order to use the quasi-pure state ap- 
proach) we shall neglect the dephasing contribution to 
73i and the feedback, through spontaneous emission, 
from state |3) to state We still allow for T3i ^ 731, 
and note the relation between the atomic dipole moment 
matrix element d\s and I^: 



r w o"l3 
t 31 



(43) 



The equations of motion for the coefficients c n ,C^ and 
C 2 are 



... 3e 

731 



(44a) 



C 3 + igi 3 



E 



i—C 2 (44b) 



C 2 = -i -fc 3 



(44c) 



where 313 = {d\z / K){hujQ / t^AL) 1 / 2 . Assuming that all 
the coefficients vary sufficiently slowly, the equation (I44bj) 
for C3 can be adiabatically integrated with the result 



77. : * C 2 {t) (45) 

731/2 - muj 731 



When this is substituted into the equation (|44c[l for C 2 , 
the decay rate f2 2 /27 3 i appears multiplying C 2 itself. 
Since this rate is, presumably, much greater than the 
transparency bandwidth (|40[) . it is consistent to assume 
that all the relevant modes are slower than it, and to 
perform a further adiabatic integration, with the result 



c n (t)e~ 



(731/2 - inuj) (^2/2731 - into) 



(46) 



Finally, this can be substituted again in (|45|) , and the 
lowest-order nonvanishing contribution in nuj kept, to 
yield the occupation probability amplitude for level 3 in 
the presence of the pulse: 



(note that this is consistent with the semiclassical treat- 
ment of [Hj], which yielded an occupation probability of 
level 3 proportional to the square of the time-derivative 
of the field amplitude envelope). 

We assume that irreversible processes, represented by 
the rate 731, take the system out of the state 3 and de- 
stroy the coherence, and we can estimate then a single- 
atom "loss" probability by 



73l|C 3 l| 2 cft' 



1673lffl 3 L 



E 



K0 2 | Cn (o)| ; 



873^31 <Ja s 2 

—q 2 — "J (ow pulso ) 



(48) 



where T = L/c is the "quantization time," Eq. (|4"3"|) has 
been used, and 5a; pu i se (the standard deviation of nu; for 
the pulse) has been defined in a natural way. Multiplying 
this by the total number of atoms, pAl, with which the 
pulse interacts, we obtain the total loss probability, 



-^loss 



5u> 



pulse 



Aw tr 



>r 2 . 



(49) 



assuming that the relevant "medium bandwidth" to be 
used in the calculations in the previous section (in partic- 
ular, in Eq. (|3~lj) ) is of the order of Acj trans , and also that 
the "effective frequency support" of the pulse, Aw pu i se , 
is of the order of a few standard deviations. 

Eq. (|4"9"|) shows that, in order to prevent the loss of 
coherence through spontaneous emission out of the level 
3, the parameter r needs to be kept very small, in which 
case, as shown in the previous section, the phase shift in 
the Schrodinger picture is necessarily very small as well. 
It may be tempting to try to look for an "optimal" pulse 
shape that, for instance, maximizes J-"o and <p (Eqs. (|32[> . 
(|33p ) while minimizing Eq. (|49|) . but that would be miss- 
ing the point. The basic meaning of Eq. (|4"9"]) is actually 
that the unitary evolution under the Hamiltonian 
assumed in the previous Section, simply does not hold 
unless the pulse's frequency spectrum is well within the 
medium's (EIT) transparency bandwidth, in which case 
r, and the maximum phase shift 0, are necessarily small. 

As an example, suppose one has a Gaussian pulse of 
the form I a ~ e -<>- z i) l a ( i n which case (<5wp U isc) 2 = 
c 2 /2a 2 , and r is given by Eq. d34j) . Then Eq. ^ be- 
comes Pioss = 2r 2 /71", and to have Pi OS s smaller than, say, 
0.1, we require r < 0.4. If we also want 1 — Tq ~ 0.1, 
Eq. (f3"2"| shows that $ cannot exceed 0.66, and then, by 
Eq. (|33|) . we have cf> < 0.26. However, this is such a small 
phase shift that the overlap of the initial state in Eq. ((T|) 
with the target state is already cos 2 (0/2) = 0.98. This 
means that one has a bigger "success probability" if one 
simply does nothing at all to the initial state. 



IV. CONCLUSIONS 



IC3I 2 



I6S13 

n 2 



nu)c n e 



(47) The results presented here are fully in agreement with 

the analysis of Shapiro and co-workers [a 0|- ^ n P ar " 
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ticular, in the "fast nonlinearity" regime, the achievable 
phase shift is very small for as long as the instantaneous 
response approximation is justified. This corresponds to 
being allowed to neglect the higher-order (in nu>) terms in 
the adiabatic expansion in Section 3, which are responsi- 
ble for the breakdown of unitary evolution as the pulse's 
bandwidth approaches the EIT transparency bandwidth. 
Note that in the formalism used in Section 3 this loss of 
unitarity is ultimately due to the disappearance of a pho- 
ton from the system (the photon is absorbed, to bring the 
atom to level |3), and then spontaneously emitted into 
some other mode) ; if this was to be described using field 
operators, restricted to only the two sets of modes "a" 
and "b" , one would have to throw in a Langevin noise 
term to preserve the commutation relations. This would 
connect to Shapiro's explanation of the reduced fidelity 
in this regime in terms of phase noise. (The connection 
is, essentially, the fluctuation-dissipation theorem.) 

A somewhat surprising result from the analysis in Sec- 
tion 2 is the realization that the seemingly arbitrarily 
large phase obtained in the Heisenberg picture does not 
necessarily translate into a large phase in the Schrodinger 
picture. This highlights an important feature of the mul- 
timode calculations. For a single mode it is certainly the 
case that any phase factor acquired by the Heisenberg- 
picture operators a(t), b(t), will also appear multiplying 
the single-photon state |11) in the Schrodinger picture. 
In the multimode case one cannot count on such a corre- 
spondence. This is already apparent from the fact that, 
in the Heisenberg picture, the magnitude of the phase 
depends on the local pulse intensity (as in Eq. [[IT])), 
whereas the Schrodinger picture treatment determines a 
single value for the phase shift, simply by projecting the 
final field state onto the initial one, as in Eq. (|2"4")) . 

It seems legitimate to say that spontaneous emission 



is ultimately responsible by the impossibility to get large 
phase shifts; in the "slow regime," as just described, by 
removing a photon from the system, and in the "fast 
regime" by populating all the initially empty pairs of 
temporal modes (with the same fj,) with equal proba- 
bility. This is consistent with many previous results that 
indicate that in order to carry a nontrivial quantum log- 
ical operation (i.e., one that can change a state into an 
orthogonal one) with an error probability P e one needs 
of the order of 1/P e photons, since 1/n is precisely the 
ratio of "spontaneous emission noise" to "signal," when 
one has n control photons [l4|. In other words, single- 
photon quantum optical gates are bound to have failure 
probabilities of the order of unity. This is plainly the 
case in schemes such as "linear optics quantum comput- 
ing" [Til ], which, however, have the advantage, over the 
schemes considered here, of not modifying the shape of 
the pulses when they succeed. 

In view of all the evidence gathered thus far, and 
the possible pitfalls of incomplete analyses, it seems rea- 
sonable to suggest that any future proposals of "single- 
photon Kerr nonlinearities" for quantum logic should at 
least include detailed studies involving: (1) clearly lo- 
cal, and physically realizable Hamiltonians; (2) localized 
wavepackets, described by quantized multimode fields, 
where at least enough modes are included in numeri- 
cal calculations to cover the whole nonlinear medium's 
bandwidth; (3) conventional fidelities computed in the 
Schrodinger picture; and (4) a realistic estimate of any 
residual losses or decoherence mechanisms. As the re- 
sults presented here indicate, however, there seems to be 
no reason to believe that such an analysis could violate 
the conclusions of Shapiro's analysis jy]. 
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